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ABSTRACT 



< 
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O . Asteroseismic data can be used to determine stellar surface gravities with preci- 

^ ' sions of < 0.05 dex by using the global seismic quantities (Av) and i^max along with 

^ , standard atmospheric data such as Tcff and metallicity. Surface gravity is also one of 

the four stellar properties to be derived by automatic analyses for 1 billion stars from 
Gaia data (workpackage GSP_Pliot). In this paper we explore seismic data from main 
sequence F, G, K stars (solar-like stars) observed by the Kepler spacecraft as a poten- 
■ tial calibration source for the methods that Gaia will use for object characterisation 

(logg). We calculate logg for some bright nearby stars for which radii and masses are 
known (e.g. from interferometry or binaries), and using their global seismic quantities 
in a grid-based method, we determine an asteroseismic logg to within 0.01 dex of the 
direct calculation, thus validating the accuracy of our method. We also find that errors 
i in adopted atmospheric parameters (mainly [Fe/H]) can, however, cause systematic 

' errors on the order of 0.02 dex. We then apply our method to a list of 40 stars to deliver 

precise values of surface gravity, i.e. uncertainties on the order of 0.02 dex, and we find 
agreement with recent literature values. Finally, we explore the typical precision that 
we expect in a sample of 400+ Kepler stars which have their global seismic quantities 
measured. We find a mean uncertainty (precision) on the order of better than 0.02 dex 
in logg over the full explored range 3.8 < logg < 4.6, with the mean value varying 
I only with stellar magnitude (0.01 — 0.02 dex). We study sources of systematic errors 

in logg and find possible biases on the order of 0.04 dex, independent of logg and 
magnitude, which accounts for errors in the Tcff and [Fe/H] measurements, as well 
as from using a different grid-based method. We conclude that Kepler stars provide a 
wealth of reliable information that can help to calibrate methods that Gaia will use, in 
particular, for source characterisation with GSP_Phot where excellent precision (small 
uncertainties) and accuracy in log g is obtained from seismic data. 

Key words: asteroseismology - stars: fundamental parameters - stars: late-type - 
surveys: Gaia - surveys: Kepler - Galaxy: fundamental parameters - 



1 INTRODUCTION 

Large-scale surveys provide a necessary homogenous set of 
data for addressing key scientific questions. Their science- 
driven objectives naturally determine the type of observa- 
tions that will be collected. However, to fully exploit the 



survey, complementary data, either of the same type but 
measured with a different instrument or of a different ob- 
servable, needs to be obtained. Combining data from several 
large-scale surveys can only result in the best exploitation 
of both types of data. 
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The ESA Gaia missiorQ is due to launch in Autumn 
2013. Its primary objective is to perform a 6-D mapping of 
the Galaxy (3 positional and 3 velocity datajf] by observing 
over 1 billion stars down to a magnitude of = 20. The 
mission will yield distances to these stars, and for about 
20/100 million stars distances with precisions of less than 
1%/10% wiU be obtained. 

Gaia will obtain its astrometry by using broad band 
G photometrj[f|. The spacecraft is also equipped with a 
spectrophotometer comprising both a blue and a red prism 
BP/RP, delivering colour information. A spectrometer will 
be used to determine the radial velocities of objects as far 
as G = 17 (typical precisions range from 1-20 kms~^), 
and for the brighter stars (G < 11) high resolution spec- 
tra (R~ll,500) will be available. 

One of the main workpackages devoted to source char- 
acterisation is GSP_Phot whose objectives are to obtain stel- 
lar properties for 1 billion single stars by using the G band 
photometry, the pa rallax tt, and the sp ectrophotometric in- 
formation BP/RP (jBailer-Jonesll2010r i. The steUar proper- 
ties that will be derived are effective temperature T^s, sur- 
face gravity logg, and metallicity [Fe/H], and also extinc- 
tion i n the astrometric G band to each of the stars. 
iLiu et al.1 (|2012l ) discuss several different methods that were 
developed to determine these parameters using Gaia data 
and we refer to this paper and references within for de- 
tails. In brief, they discuss the reliability of determining 
the four parameters by using simulations, and in particular, 
they conclude that they expect typical precisions in logg 
on the order of 0.1 - 0.2 dex for main sequence late-type 
stars, with mean absolute residuals (true value minus in- 
ferred value from simulations) no less tha n 0.1 dex f o r star s 
of all magnitudes, see Figure 14 and 15 of lLiu et al.l (|2012l) . 
We note that the stellar properties derived by GSP_Phot will 
be used as initial input parameters for the workpackage de- 
voted to detail ed spectroscopic analysis of the brighter tar- 
gets GSP_Spec (iRecio-Blanco et aIll2006l:lBiiaoui et al.l l 2010l : 



iKordopatis et al.ll201ll : IWorlev et al]l2012h using the Radial 
Velocity Spectrometer data (|Katj|2005l v Spectroscopic de- 
terminations of logg, Tcff and [Fe/H] in general can have 
large correlation errors, and if logg is well constrained, Tcff, 
[Fe/H] and chemical abundances can be derived much more 
precisely. 

The differe nt algorithms for GSP_Phot discussed by 
iLiu et al.1 l|2012h used to determine the stellar properties in 
an automatic way have naturally been tested on synthetic 
data. However, to ensure the validity of the stellar proper- 
ties, a set of 40 bright benchmark stars have been compiled 
and work is still currently underway to derive stellar proper- 
ties for all of these in the most precise, homogenous manner 
(e.g. Heiter, et al. in prep.). Unfortunately, these benchmark 
stars will be too bright for Gaia, and so a list of about 500 
primary reference stars has also been compiled. The idea is 
to use precise ground-based data and the most up-to-date 
models (known from working with the benchmark stars) to 
determine their stellar properties as accurately (correct) and 



^ http:/ /sci.osa.int/scioncc-c/www/arca/indcx.cfm?faroaid=26 

^ For stars fainter than V~17, the radial velocities will not be 

available. 
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Figure 1. High SNR power /frequeney spectrum of a Kepler 
solar-like star KIC 6603624. Some individual frequencies v are 
denoted with their degre e I and radial order n in parenthesis 
llAppourchaux et al.ll201^ . The reference radial order no usually 
corresponds to orders similar to 20 for these stars. We also show 
the individual large and small frequency separations, Av(Z, n) and 
5v{l, n) respectively, for three cases, and the approximate fre- 
quency corresponding to the maximum power in the spectrum, 
i^max- For much lower SNR spectra, the mean value (Ai/) of the 
large frequency separations and I'max can usually be determined 
even if individual frequencies v can not be resolved. 



precisely (small uncertainties) as possible. These primary 
reference stars will be observed by Gaia and thus will serve 
as a set of calibration stars. A third list of secondary ref- 
erence stars has also been compiled. These consist of about 
5000 fainter targets. 

In the last decade or so, much progress in the field 
of asteroseismology has been made, especially for stars ex- 
hibiting Sun-like oscillations. These stars have deep outer 
convective envelopes where stochastic turbulence gives rise 
to a broad s pectrum of resonant osc il lation modes (e.g. 
UlrichI Il97d iLeibacher fc SteinI Il97ll. iBrown fc Gillilandl 



1994 ISalabert et al.1 |2002| . iBouchv fc Carrieij |2002| ). The 



power spectrum of such stars can be characterised by some 
global seismic quantities; (Au), i^max, and (Sv). The quantity 
(Ai^) is the mean value of the large frequency separations 
A.i'i,n ~ i^i,n — vi,n-i whcrc fi^ri is a resouaut oscillation 
frequency with degree I and radial order n, fmax is the fre- 
quency corresponding to the maximum of the bell-shaped 
frequency spectrum, and (Su) is the mean value of the small 
frequency separations Sui^n = J^i,n — i'i+2,n-i- Figure [1] shows 
the power (frequency) spectrum of a solar-like star, observed 
by the if ep/er spacecraft, depicting these quantities. 

Even when individual frequencies can not be determined 
from the frequency spectra both (Av) and z^max can still 
be extracted quite robustely. Many methods have been de- 
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mic quantities have been shown to scale with stellar param- 
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The G photometric scale is similar to V. 



Silva Aguirre et al.l (|2012l ). 

seismic quantities with the observed ones over a large grid of 
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stellar models, very precise determinations of \ogg (< 0.05 
dex) and mean density (< 2%) can be obtained for main 
sequence F, G, K star s (iBruntt et al.' '201C ; iMetcalfe et al 
20ld: iGai et all I2OIII : iMathur et al., ,2013 : ICreevev et"ar 



Table 1. Summary of tlie frequently used notation in tliis study. 



2012al ). 



Of particular interest for Gaia is the Keple 70 field of 
view, ~100 square-degrees, centered on galactic coordinates 
76.32°, +13.5°. Kepler is a NA SA mission dedicated to char- 
acterising planet-habitability jBorucki et al.l [2OI0I I. It ob- 
tains photometric data of approximately 150,000 stars with 
a typical cadence of 30 minutes. However, a subset of stars 
(less than 1000 every month) acquire high-cadence data with 
a point every 1 minute. This is sufficient t o detect and 
charac teris e Sun-like oscillations in many stars. IVerner et al] 
(|201ll ) and lChaplin et al.1 (|201ll ') recently showed the detec- 
tions of the global seismic quantities for a sample of 600 F, 
G, K main sequence and subgiant (V /TV) stars with typical 
magnitudes 7 < V < 12, while both CoRoT and Kepler have 
both shown their capabilities of detecti ng these same seismic 
quantities in 1000s of red giants, e.g. [Hekker et al.l (l2009f): 
[l<;allingcr ct al.' ('2010''): 'Bedding ot all (I2OIOI'): iBaudin et all 
(2011); Mosscr ct al . (2012); St cUo ct al, ( 20131 ). 

With the detection of the global seismic quantities in 
hundreds of main sequence stars (and 1000s of giants), the 
Kepler field is very promising for helping to calibrate Gaia 
GSP_Phot methods. In particular, the global seismic quan- 
tities deliver one of the four prop erties to be extr acted by 
automatic analysis, namely log q. [Gai et al.l (|201ll ) studied 
the distribution of errors for a sample of simulated stars us- 
ing seismic data and a grid-based method based on stellar 
evolution models. They concluded that \ogg derived from 
seismic properties {"seismic logg") is almost fully indepen- 
dent of the input physics in the stellar evolution m odels 
that are used. More recently, iMorel fc Miglid (|2012l ) com- 
pared classical determinations of log g to those derived alone 
from a scaling relation (see Eq.[2]), and concluded that the 
mean differences between the various methods used is ~0.05 
dex, thus supporting the validity of a seismic determination 
of logg. While some studies have focussed on comparing 
seismic r adii or masses with alternative determinations, for 
example, iBruntt et al.l (|20ld ). no study has been done fo- 
cussing on both the accuracy and precision of a seismic log g 
using one or more grid-based methods for stars with inde- 
pendently measured radii and masses. The accuracy and 
precision in logg for these bright stars has also not been 
tested while considering precisions in data such as those ob- 
tained by Kepler. Such a study could validate the use of 
seismic data as a calibration source for Gaia. For the rest of 
this work the term precision refers exclusively to the derived 
uncertainty in logg, while accuracy refers to how true the 
value is. 

With these issues in mind, the objectives of this paper 
are to (i) test the accuracy and precision of a seismic logg 
from a grid-based method using bright nearby targets for 
which radii and masses have been measured (Sect. 3), (ii) 
determine logg for an extended list of stars whose global 
seismic properties and atmospheric parameters are available 
in the literature using the validated method (Sect. 4), and 
(iii) study the distribution of log g and their uncertainties of 
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Notation 


Definition 
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gravity 


logs 


logarithm of g 


TcS 


effective temperature 


(Au) 


mean large frequency separation 


t'max 


frequency of maximum power 


[Fc/H] 


mctallicity 


R 


stellar radius 


M 


stellar mass 


cr 


uncertainty in log g 


s 


systematic error in log g 


Ac 


extinction in the Gaia G band 


r 


stellar magnitude in the SDDS system 



over 400 F, G, K V/IV Kepler stars derived by a grid-based 
method while concluding on realistic uncertainties (preci- 
sions) and possibles sources of systematic errors for this po- 
tential sample of Gaia calibration stars (Sect. 5). In this 
work, our analysis is restricted to stars with logg > 3.75 
primarily due to the very limited sample of stars for which 
we can test the accuracy of and validate our method. For 
ease of reading. Table [T] summarizes the frequently used no- 
tation in this work. We begin in Sect. 2 by summarising the 
different methods available for determining logg. 



2 DIRECT METHODS TO DETERMINE log g 

In this section we summarise the various methods that are 
used to determine the surface gravity (or the logarithm of 
this logg) of a star. Comparing each of these methods di- 
rectly would be the ideal approach for unveiling shortcom- 
ings in our models (systematic errors) and reducing uncer- 
tainties by decoupling stellar parameters. 



2.1 Derivation of logg from independent 
determinations of mass and radius 

The most direct method of determining log g involves mea- 
suring the mass M and radius J? of a star in an independent 
manner. Surface gravity g is calculated using Newton's Law 
of Gravitation: g = GM/R^ where G is the gravitational 
constant. 



2.1.1 Mass and radius from eclipsing binary systems 

For detached eclipsing spectroscopic binaries, both M and 
R can be directly measured by combining p hotometric and 
radia l velocity time series jRibas ct al... 20051 : ICreevev et al.l 
l2005l . I2OIII : iHelminiak et al.1 l2012l l~The orbital solution is 
sensitive to the mass ratio and the individual Msini of 
both components, where i is the inclination. The photomet- 
ric time series displays eclipses (when the orbital plane has 
a high enough inclination) that are sensitive to i and the 
relative R. Once i is derived, then the individual M axe 
solved. Kepler's Law relates the orbital period of the system 
n, the system's M, and the separation of the components, 
n is known from either eclipse timings or observing a full 
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radial velocity orbit. Once the individual M are known then 
n scales the system (providing the separation) and thus the 
individual R, and g can be then calculated. 



2.1.2 Mass and radius from interferometry and global 
seismic quantities 

R is measured by combining the angular diameter 6 from in- 
terferometry with the distance to the star. The distance (or 
its inverse the parallax) has been made available using data 
from the Hipparcos sa t ellite for stars with V <~ 8 mag 
llPerrvman et al.l 1 19971 : Ivan LeeuwenI l2007l : iKervella et all 
Indirect methods also exist for determining the an- 
gular diameter of a star, such as combining Tcs with m ea- 
surements of bolometric flux (|Silva Aguirre et al.ll2012|). or 
from calibrated relations using photometry ( Kervella et al.l 
l2004h . Using R and {Av) a model-independent mass deter- 
mination can be obtained by using the asteroseismic relation 
which links mean density (p) and {Au) 



IJL = ^[m/Mq)/{r/RqY 



(1) 



where (Ai')(^ is the so lar value (e.g. iKieldsen fc Bedding! 
1 19951 : iHuber et al.ll201ll ) . 



2.1.3 Mass and radius from interferometry and/or high 
precision seismic data 

When high signal-to-noise ratio seismic data are available, 
individual oscillation frequencies (see Fig. [T| can be used 
to do detailed m o delling, and hence d etermine M (e.g . 



Dogan et al.1 12OI0I: iBrandao et al.l I2OIII: iBigot et all I2OIII : 



Metcalfe et al.l 



2012al : iMiglio et al.H2oTa) . When combined 
with an independently measured R, the uncertain t ies in 
mass a (M) can reduce to < 3% (jCreevev et al.l l2007l : 



iBazot et al]|201ll : iHuber et aLll2012l ). If an independently 
determined R is not available, then stellar modelling also 
yields R to between 2 - 5% and M with less precision. How- 
ever, this method depends on the physics in the interior stel- 
lar models unlike the methods mentioned above, and using 
different input physics may result in different values of mass. 
Typical uncertainties/accuracies in M does not usually ex- 
ceed about 5% for bright targets when more constraints are 
available, and this translates to less than a 0.02 dex error in 
logg for stars that we consider in this work. As g is sensi- 
tive to seismic data, then R and M are correlated and log g 
can de determined with a precision of ~ 0.02 dex with a 
very slight dependen ce on the physics in the models, e.g. 



very sUgnt dependen ce 
iMetcalfe et ail (|2010l ). 



2.2 Spectroscopic determinations of log 3 

The surface gravity of a star is usually derived 
from an atmospheric analy si s wit h spectroscopi c data 
fe.g. iThevenin fc Jasniewicj Il992l: iBruntt et al.1 I2OI0I : 
iLebzelter et al.ll2012l : ISousa et al.ll2012l ). There are two usual 
approaches for deriving atmospheric parameters (logg, Tcff, 
and [Fe/H]). The first approach is based on directly com- 
paring a library of synthetic spectra with the observed one, 
usually in the form of a best-fitting approach. A shortcom- 
ing of this method is that combinations of parameters can 



produce similar synthetic spectra so that many correlations 
between the derived parameters exist. The more classical 
method for determining atmospheric parameters relies on 
measuring the equivalent widths of iron lines (or other chem- 
ical species). This method assumes local thermodynamic 
equilibrium (LTE) and requires model atmospheres. Once 
Tefi is determined (by requiring that the final line abun- 
dance is independent of the excitation potential or for stars 
with rcff> 5000 K, measuring the Balmer H line profiles), 
then g is the only parameter controlling the ionisation bal- 
ance of a chemical element in the photospheric layers, which 
acts on the recombination frequency of electrons and ions. 
g is then determined by requiring a balance between differ- 
ent ionized lines e.g. Fe-I and Fe-II. Spectroscopically de- 
termined logg can have large systematic errors especially 
for more metal-poor stars where NLTE effects must be 
taken into account. This can change loggr by '^0.5 dex, e.g. 
IThevenin fc Jasniewicj (|l992l) . 



2.3 \ogg derived from Hipparcos data 

An alternative method for determining log g relies on know- 
ing the distance d (in pc) to the star from astrometry, its 
bolometric flux (combining these gives the luminosity of a 
star) and Tcff. Then substituting M and logg for R in Ste- 
fan's Law, one obtains the following logg = —10.5072 -|- 
logM-|-0.4(Vo-|-BCv)-21ogd-|-41ogreff, where Vo is the de- 
reddened V magnitude and BCy the bolometric correction 
in V, and M (in M0) a fixed value th at is estimat ed from 
evolution models, e.g. IThevenin et al.l | 2001 : Barbuv et al.l 
I2OO3I . This Hipparcos log g is often used as a fixed parameter 
for abundance analyses of stars. Typical uncertainties are no 
less than 0.08 dex where especially the error in M is large, 
and errors from d and TcS are not insignificant. We note 
that Gaia will deliver unmistakeably accurate distances for 
much fainter stars and these will provide a much improved 
log g using this method for individual stars. 



2.4 log g from evolutionary tracks 

2.4.1 Classical constraints m the H-R diagram 

When H-R diagram constraints are available (Toff, L, metal- 
licity) then stellar evolution tracks can be used to pro- 
vide estimates of some of the other parameters of the star, 
e.g. mass, radius, and age. While correlations exist between 
many parameters, e.g. mass and age, these correlations also 
allow us to derive certain information with better precision, 
e.g. mass and radius gives g. Exploring a range of models 
that pass through the erro r box thus allows us t o limit the 
possible range in logp (e.g. ICreevev et al.|[2012bl ). 



2.4.2 Grid-based asteroseismic logg 

Apart from performing detailed modelling using asteroseis- 
mic data (Sect. [27l.3|) . one can rely on grids of stellar models 
to estimate stellar properties such as mass and radius with 
precisions on the order of 2-12%. However, because astero- 
seismic data are extremely sensitive to the ratio of these two 
parameters, then very precise determinations of \ogg and 
(p) can be obtained in a n almost model-independent man- 
ner (e.g. iGai et al.|[2oTll . in Sect. 15.31 we address this issue). 
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Such a grid-based asteroseismic \ogg can be obtained in the 
following manner: a large grid of stellar models that spans 
a wide range of mass, age, and metallicity is constructed. 
Each model in the grid has a corresponding set of theoret- 
ical observables, such as logg, Tefi, individual frequencies 
1^71,1, (Ai/) and z^max- A scaling relation is used to obtain 

i^max ■ 

(M/Mo) 



, = (2) 

where t^max.a is the solar value (e.g. IVandakurovl 1 19681 : 
lTassouJ[li80l : lKieldsen fc Beddinjl995h . The quantity (Ai/) 
can be obtained either from Eq. [1] or calculated directly 
from the oscillation frequencies derived from the structure 
model. Differences, however, on the order of 2% in (A^) 
may be found by adoptin g one or the other method (e.g., 
IWhite et al]|201ll . see also lMosser et al.1 120131 '). A set of in- 
put observed data, e.g. {(Ai/), Tetr, and [Fe/H]}, is com- 
pared to the theoretical one, and the models that give the 
best match to the data are selected and the value of \ogg 
and its uncertainty is derived using these best selected mod- 
els. Such a precise value for this seismic logp comes, in fact, 
from the very close relation between i^max and the cut-off 
frequency, and recent work has made progr ess on under- 
standing this relation (|Belkacem et al.ll201ll ). Several seis- 
mic grid-based approaches for determining stellar proper- 
ties have been discussed and applied in the recent litera- 



ture dStello et al.ll2009l:lQuiriori et al.ll2010l:lBasu et a,l.ll201C 



Gai et al.ll201ll: ICreevev et al.ll2012a[ : lMetcalfe et al.ll2012~ 
Silva Aguirre et al.ll201 j T and for the rest of this paper we 
use the term seismic logg to refer specifically to the grid- 
based method for determining log p. 



3 COMPARISON OF DIRECT AND 
SEISMICALLY DETERMINED log g. 

3.1 Observations and direct determination of logi; 

In order to test the reliability of an asteroseismically deter- 
mined logg, the most correct method is to compare it to 
log g derived from mass and radius measurements of stars in 
eclipsing binaries (Sect. 2.1.1), apart from the Sun. Unfortu- 
nately the number of stars whose masses and radii are known 
from binaries, where seismic data are also available, is quite 
limited. For this sample of stars we have a Cen A and B, 
and Procyon. Following this, we rely on the combination of 
asteroseismology and interferometry to determine logg, and 
use the scaling relation which links density to the observed 
properties of (Az/) and R to provide an independent M mea- 
surement (Sect. 2.1.2). However, since this scaling relation is 
used explicitly in the grid-based method, we have opted to 
omit stars where this method provides the mass, except for 
the solar twin 18 Sco, which was included because of its sim- 
ilarity to the Sun. To complete the list of well-characterised 
stars we have then chosen some targets for which an interfer- 
ometric R has been measured and detailed seismic modelling 
has been conducted to determine the star's M by several au- 
thors (Sect. 2.1.3). The stars which fall into this group are 
HD 49933 and /? Hydri. The seven stars are listed in order 
of logg in Table [2] along with (Ai/), i/^ax. Toff, [Fe/H], M, 
and R. When several literature values are available these 



are also listed. The final column in the table gives the direct 
value of log g as derived from M and R. For HD 49933 and 
P Hydri we adopt the weighted mean values of log g which 
are 4.214 and 3.958 dex respectively (see Table [2]| and these 
are summarized in column 2 of Table [3] For the rest of the 
paper we refer to these determinations of log g as the direct 
determinations. 



3.2 Seismic method to determine log g 

We use the grid-based method, RadEx lO, to deter- 
mine an asteroseismic value of logg (|Creevev et al.l 
l2012ah . The grid was constructed using the A STEC 
stellar evolution code (|Christensen-Dalsgaardl |2008| ) with 
the following configu ration: the EFF equation of state 
(|Eggleton et al.l Il973h without Coulom b corrections, the 
OPAL opacities ( Iglesias fc Roger3ll996l ) supplemented by 
Kurucz opacities a t low temperatures, solar mixture from 



Grevesse fc Noelsl (Il993h and nuclear reaction rates from 
Bahcall fc Pinsonneault (| 19921 ). Convection in the outer con- 



vective envelope is de scribed by the mixing-length theory of 
iBohm-Vitensd (|l958h and this is characterised by a variable 
parameter omlt (where I = otyiurHp, I is the mixing-length 
and Hp the pressure scale height). When a convective core 
exists, there is an overshoot layer which is also characterised 
by a convecctive core overshoot parameter Qov and this is 
set to 0.25 (an average of values used recently in the lit- 
erature). We also ignore diffusion effects, although we note 
that for accurate masses and ages this needs to be taken into 
account. 

The grid considers models with masses M from 0.75 - 
2.00 M© and ages t from ZAMS to subgiant. The initial 
metallicity Zi spans 0.003 - 0.030 in steps of ^ 0.003, while 
initial hydrogen abundance Xi is set to 0.71: this corresponds 
to an initial He abundance Vi = 0.260 — 0.287. The mixing 
length parameter omlt = 2.0 is used, which was obtained 
by calibration with solar data. We note that to derive other 
stel lar properties, i t is very important to a l low q mlt to vary 
e.g. ICreevev et al] (|2012bii : iBonaca et all (|2012D . 

To obtain the grid-based model stellar properties, e.g. 
logg, M, and t, we perturb the set of input observations 
by scaling its input error with a number drawn randomly 
from a normal distribution and add this to the input (ob- 
served) value. We compare the perturbed observations to the 
model ones and select the model that matches best. This is 
repeated 1,000 times to yield model distributions of best- 
matching stellar parameters. In this method, (Ai/) is calcu- 
lated using Eq.[T]and the input observations consist primar- 
ily of (Al/), z^max, Tcff, aud [Fe/H], although other inputs are 
possible, for example, L or R. The model distributions are 
then fitted to a Gaussian distribution and the fitted stellar 
property and its uncertainty {a) are defined as the central 
value and the standard deviation a of the Gaussian fit. In 
this work we consider just the derived value of logg. 



3.3 Analysis approach 

We determine a seismic logg for the reference stars using 
the method explained above. We consider different sets 
of input data in order to test the effect of the different 
observations on the accuracy (and precision) of a seismic 
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Table 2. Properties of the reference stars 



Star 


{Au) 


^max 




[Fe/H] 


R 


M 


logg 




(/i,Hz) 


(mHz) 


(K) 


(dex) 


(Rg) 


(Mq) 


(dex) 


«CenB 


161.5 ± 0.11^" 


4.0^" 


5316 ± 28" 


+0.25 ± 0.04" 


0.863 ± 0.005" 


0.934 ± 0.006" 


4.538 ± 0.008 












0.863 ± 0.003" 






18 See 


134.4 ± 0.3^" 


3.l2" 


5813 ± 21^" 


0.04 ± 0.01^" 


1.010 ± 0.009^" 


1.02 ± 0.03^" 


4.438 ± 0.005 


Sun 


134.9 ± 0.1^" 


3.05="' 


5778 ± 20^= 


0.00 ± 0.01 






4.437 ± 0.0023<* 


aCcnA 


105.6*" 


2.3*" 


5847 ± 27" 


+0.24 ± 0.03" 


1.224 ± 0.003" 


1.105 ± 0.007" 


4.307 ± 0.005 



2 44!, 

HD49933 85.66 ± 0.18^" 1.8^" 6500 + 75^'' -0.35 + 0.10^'' 1.49^'' 1.325^'' 4.23 + 0.025'' 

1.657 + 0.028^'= 

6640 + loos'* -0.38^'' 1.42 + 0.045<* 1.20 + 0.08^'' 4.212 + 0.039 

1.39 + 0.045*= + 0.035^ 4.201 + 0.027 

1.44 + 0.04^= 1.20 + 0.035= 4 200 + 0.027 

Procyon 55.5 + 0.5'^" l.O'^'' 6530 + 90*''= -0.05 + 0.03^^* 2.067 + 0.028^= 1.497 + 0.037''-'' 3.982 + 0.016 

/3Hydri 57.24 + 0.16^" 1.0^" 5872 + 44^* -0.10 + 0.07'^= 1.814 + 0.017'^'' 1.07 + 0.03'''' 3.950 + 0.015 

1.04^'' 3.938 + 0.015 

1.082''= 3.955 + 0.015 

l-10^o'.07 ^■'^ 3.962 + 0.030 

5964 + 70''9 -0.03 + +0.07''9 1.69 + 0.05''s 1.17 + 0.057f 4.02 + 0.0479 



Re ferenc es: HKieldsen et al.l ll2005h . ^ IPorto de Mello et al.l l|2008l ). ^ feervella et al.l ll2003l) . -^'IPourbaix et al.l ||2002|'I. HBmot et al.l 
||2006|'I. ^"iBazot et al.l ll2qill'), 3" Taking the average of Table 3 from l^utain fc ProehlichI lll992l '). ^i jKieldsen fc Beddintd jl99'5li ^ 
3t|Qj.gygggg Sauvall lll998l'l. 3''lMm=1 . 98919 e30 kg, lR0=6.9599e8 km, G = 6.67384ell kg-i s-2(NIST database), 
*°lBouchv fc Carried (120021). *lQuirion et all boiol ). 5" Using the Z = modes with Height/Noise>l from Table 1 of iBenomar et al.1 
20091). S iKaUinger et al.l 1I2OIOI) log g is the mean and sta ndard deviation of the val ues cited from their Ta ble 1. S'iGru berbauer et al.l 



1 I2OO9I), S^Bigot et al.l ll201l[). ^'ICreevev fc Bazoli 1I2OIII). ^"lEggenberger et al.l ll2004l ), ^'jMartic et al.l f2004l ). '^ 'iFuhrmann et al.l 



S'lAllende Prieto et al 1I2OO2I). ^^Kervella et al.l ll2004^^J^Girard et alT i 2000l). ''"iBedding et"aiTi200'^ . ' ^iNorth et al.l (120071)!" 
^ teruntt et al.l ll2010l )~' 1Brandao et al.l 1 I2OIII) . ' HPogan et al.l | |2010^ ■ '^ ^F^nandes fc Monteird | |2003^ . '^ ^da Silva et al.l l l20od ) 



(51) {(Ai.), i.n.ax, Teff, [Fe/H]}, 

(52) {(Al.), i^^^ax, Teff}, and 

(53) {(Al.), Teff}, 

For the potential sample of Gala calibration stars, 
[Fe/H] will not always be available, and in some cases, i^max 
is difficult to determine for very low S/N detections. 

The observational errors in our sample are very small 
due to the brightness (proximity) of the star, so we also 
derive an asteroseismic log g while consid ering observational 
errors that we expect for i^ep/er stars (see lVerner et allboill 
and Figure [SJ. We consider three types of observational er- 
rors while repeating the exercise: 

(El) the true measurement errors from the literature, 
(E2) typically "good" errors expected for these stars, i.e. 
o-((Ai/» = 1.1 /iHz, a(f„,ax) = 5%, o-(reff) = 70 K, and 
a([Fe/H]) = 0.08 dex (see Sect. [53]), and 
(E3) "not-so-good" measurement errors, primarily consider- 
ing the fainter targets (V~ll,12), i.e. (t{{/^v)) = 2.0 /iHz, 
(T(//,nax) = 8%, cr(reff) = 110 K, and cr([Fe/H]) = 0.12 dex. 



3.4 Seismic versus direct log 17 

In Figure [5] we compare the grid-based seismic log g with 
the direct logg for the seven stars. Each star is represented 
by a point on the abscissa, and the y-axis shows (seismic - 
direct) value of logg. There are three panels which repre- 
sent the results using the three sets of input data. We also 
show for each star in each panel three results; in the bot- 



tom left corner these are marked by 'El', 'E2', and 'E3', 
and they represent the results using the different errors in 
the observations. The black dotted lines represent seismic 
minus direct \ogg = 0, and the grey dotted lines indicate 
±0.01 dex. 

Figure [5] shows that for all observational sets and errors 
logg is generally estimated to within 0.01 dex in accuracy 
and with a precision of 0.015 dex. This result clearly shows 
the validity of the global seismic quantities and atmospheric 
parameters for providing accurate values of log g. It can also 
be noted that the precision in log g typically decreases as (1) 
the observational errors increase (from E2 - E3), and (2) the 
information content decreases (SI - S2 - S3, for example). 

One noticeable result from Figure [2] is the systematic 
offset in the derivation of log g for HD 49933 when we use 
[Fe/H] as input (SI). This could be due to an incorrect 
[Fe/H] or T^s, an error in the adopted direct log <; or a short- 
com ing of the grid of mo dels. This star is known to be active 
and lMosser et all (|2005l ) found clea r evidence of s pot s igna- 
tures in line bisectors . More recentlv lGarcfa et all (I2OIOI ') and 
ISalabert et al.l (|201ll ) found frequency and amplitude vari- 
ations similar to those in the Sun, ev idence of the presence 
of a stellar cycle. I Garcia et all (|2010l ) derived an S-index of 
about 0.3 corresponding to a very active star, even though 
no detection of a magnetic field has yet been confirmed (P . 
Petit, private communication 2012). iFabbian et al.l (|2012l ') 
showed that the magnetic effects in stellar models affect the 
determination of the metallicity and this, in turn, will affect 
the determination of other stellar parameters such as mass. 
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Figure 2. Accuracy of method. Seismic-minus-direct logg for the seven sample stars while considering different subsets of input 
observations (different panels, SI, S2, S3) and different observational errors (El, E2, E3). See Sect. 1X51 for details. 



In their work, they considered large magnetic field strengths 
of several hundreds of Gauss. 

Since no evidence of such a large magnetic field has yet 
been found for HD 49933 we propose another explanation. 
We consider the effect of the confirmed presence of spots 
on the effective temperature. By considering a spot area of 
about 20 percent of the stellar surface we found that the real 
effective temperature (that of the non spotted surface of the 
star) should increase by about 300 K with respect to the 
non-spotted star. Adding 300 K to Teff results in a seismic 
logp that increases by 0.011 dex for SI (considering errors 
El and E2), which makes the new value consistent with the 
direct value within error bars. The increase of 0.011 dex 
corresponds to a relative increase of 3.7a and 1.2a for both 
error types El and E2. If we add this same 300 K for cases S2 
and S3, we also find an increase in logg, but smaller (0.009 
and 0.006 dex), corresponding to a relative increase of 0.6a 
and 0.4a, respectively. 



3.5 Systematic errors in observations 

For all of the calculations we have assumed that the input 
observations are correct (accurate). While this is certainly 



Table 3. logg derived by RadExlO for the reference stars using 
the true measurement errors. A3 = logg — loggdircct- 



Star 


logs 


log ^direct 


A3 


A3 




(dex) 


(dex) 


(dex) 




a Cen B 


4.527 ± 0.004 


4.538 


-0.011 


-2.8 


18 Sco 


4.441 ± 0.004 


4.438 


0.003 


0.8 


Sun 


4.438 ± 0.001 


4.437 


0.001 


1.0 


a Cen A 


4.312 ± 0.004 


4.307 


0.005 


1.3 


HD 49933 


4.195 ± 0.007 


4.214 


0.019 


2.7 


Procyon 


3.981 ± 0.006 


3.982 


-0.001 


-0.2 


/3 Hydri 


3.957 ± 0.010 


3.958 


0.001 


0.1 



more true for brighter nearby targets where high SNR data 
can be obtained, the same cannot be said for fainter stars. 
In particular with spectroscopic data, the determination of 
Tcfi and [Fe/H] are correlated and depend on the analysis 
methods used and the different model atm ospheres (see e.g. 
ICreevev et al ]|2012al : lLebzelter et al.ll2012l ). Additionally for 
many stars a photometric temperature may be the only 
available one and while these estim ates are very good, sys- 
tematic errors are still unavoidable (jCasagrande et aLll2010l : 
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Figure 3. Seismic-minus- direct values of logg for /3 Hydri when 
the observations witli systematic errors included for cases 2 and 
3 are analysed by the seismic method. See Sect. 13.51 



true observational errors. We highlight the excellent agree- 
ment between our seismically determined parameters, and 
those obtained by direct mass and radius estimates, log g is 
matched to within ~0.01 dex, and 0.02 dex for HD 49933. 
We must also comment on the very small uncertainties given 
in Table O these results were obtained using the true (very 
small) observational errors given in Table [2] and typically 
one would not expect to obtain such small errors for fainter 
stars. As can be seen, the results using the relaxed observa- 
tional errors give more reasonable parameter uncertainties, 
i.e. 0.010 - 0.015 dex, see Fig. (2] while still matching the 
direct logg to within la. 



iBovaiian et~aLll2012h . in particular due to unknown redden- 
ning. Larger photometric errors also lead to larger errors 
on the temperatures. This is not only a problem for fainter 
stars. For example, for j3 Hydri we found two determinations 
of Teff — an interferometric one and a spectroscopic one. For 
the spectroscopic Teff, the corresponding fitted metallicity 
from the atmospheric analysis will correlate with it. 

To study the effect of systematic errors in the obser- 
vations, we repeated our analysis for /? Hydri while using 
three sets of input data that change only in T^s and [Fe/H]. 
The first set (1) u ses the interferometric value of Teff from 
iNorth et all (|2007l ) and their [Fe/H] (5872, -0.10), which we 
consider as the most correct, the second set (2) uses (5964, 
-0.10), and the third set (3) uses (5964, -0.03) as given by 
Ida Silva et"al] (120061 ). The results for SI and S2 are shown 
in Figure O The lower panel shows that for a systematic 
error in both Teff and [Fe/H] the accuracy decreases. The 
top panel shows that when we only consider the T^s infor- 
mation we get a smaller increase in the offset than when we 
consider both T^s and [Fe/H]. One way of interpreting this 
result is by considering that in S2 there is much more weight 
assigned to the seismic data than the atmospheric data, and 
so an incorrect atmospheric parameter should not influence 
the final result as much as in case SI where the atmospheric 
parameters have more weight. In the latter case, an incorrect 
Teff with the correct [Fe/H] will necessarily shift the mass 
either up or down (in terms of the H-R diagram), and result 
in a more displaced log;;. However, in both cases, we see 
that the offset does not exceed 0.02 dex. For determining a 
seismic radius and mass, however, a systematic error in the 
input observations has a much more profound effect on the 
offset. In this case biases on the order of up to 6% and 20% 
in ra dius and mass can be obtained (jCreevev fc TheveninI 
l2012ll . It must also be noted that a systematic error in the 
atmospheric parameters is going to have a much larger neg- 
ative effect when we use only the global seismic quantities 
instead of performing a detailed seismic analysis with indi- 
vidual frequencies, where the latter have much more weight 
in the fitting process. 

3.6 Seismic determination of \ogg from the global 
seismic quantities for the reference stars 

We summarize logg for the sample stars in Table [3] derived 
by RadExlO using (Ai/), i^max, T^b, and [Fe/H], and the 



4 DETERMINATION OF log g FOR AN 
EXTENDED LIST OF STARS 

We apply our grid-based method to an extended list of 
stars with measured global seismic quantities and atmo- 
spheric parameters. Table [4] lists the star name along with 
the other measured parameters that are used as the input for 
our method. The first part of the table comprises primarily 
bright stars whose oscillation properties have been measured 
either from ground-based or spaced-based instrumentation 
(see references given in the table). For most of these stars 
no errors are cited for {Au) and i/max. The second part of 
the table lists a set of 22 solar-t ype stars obse r ved b y the 
Kepler spacecraft and studied in iMathur et all (|2012l) . We 
have taken the seismic and atmospheric data directly from 
this paper. To conduct a homogenous analysis of these stars 
we adopted a 1.1 /^Hz error on (Ai/) for all of the stars, 
and a 5% error on !^max, typical of w hat has been found for 
the large sample of Kepler stars fsee lHuber et al.ll201ll and 
Sect. [5T|) . 

Table [5] lists the derived value of log g and 2a uncertain- 
ties (to allow for round-up error) for each of the stars using 
RadExlO. We show two values of \ogg; the first is obtained 
by using the four input constraints { (Aj^),i^niax, Teff , [Fe/H]} 
and the second logpno[pe/H] is obtained by omitting [Fe/H] 
from the analysis. 

Figure |4] top panel compares the derived values of \ogg 
for the Kepler stars with those determi ned using the individ - 
ual oscillation frequencies, as given bv lMathur et al.l (|2012h . 
with our la error bars overplotted. We see that the grid- 
based method provides \ogg consistent with those derived 
from a detailed asteroseismic analysis, although a very small 
trend can be seen. For some of the stars they obtain a fitted 
initial He abundance significantly below the accepted pri- 
mordial value, suggesting that the corresponding fitted mass 
and radius may be slightly biased (grey squares) . If we omit 
these stars, then we fit a slope of -0.03 ± 0.02 to the differ - 
enc e between thei r and o ur log g values. I White et al.] (|201ll ) 
and iMosser et all (j2013l 'l point out that a discrepancy may 
exist between the two different theoretical values of (Au) 
(using scaling relations or from individual frequencies), and 
the two different approaches could also be responsible for 
this trend. We see, however, that our values agree generally 
to within 0.01 dex. ISilva Aguirre et al.l (|2012l 'l have analysed 
6 of these stars and the lower panel of Fig. |4] shows a com- 
parison between their log g values with ours. Fitting the dif- 
ferences between our results, we obtain a slope of 0.001 ± 
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Table 4. Measured seismic and atmospheric properties for an extended list of Sun-like oscillators 



Star Name 


(A,.) 






[Fc/H] 






(mHz) 




(K) 


(dex) 


HD 


10700 


169. 0*^ 


45QQa 


c:qqq -L A7a 
0000 ^ '-t I 


1 07 


HD 


17051 


-LZU.U 


97nn^ 

Z i uu 


fiOSO -1- RO 
OUoU ^ ou 


n 1 -u n 07 

U.IO ^ U.Ul 


HD 


23249 




/ uu 


AQQfi -U K7d 

^yoo ^ 1 


n 1 -1- 07 

u. 10 ^ U.U i 


HD 


49385 


uO.o 


1 ni '^^ 


uuyo ^ uu 


in OQ -1- 0^^^ 
T^u.uy ^ u.uo 


HD 


52265 


yo.o m u.i 


1 snn^ 

loUU 


«i 00 -1- fiO^ 
D-LUU It OU 


1 Q -1- Of^^ 

u.iy m u.uo 


HD 


61421 


•JO.KJ 


1 nnn 

-LUUU 


RAQA -\- 


01 -1- 07 

U.Ul ^ U.U ( 


HD 


63077 


Q7 
y / .u 


ZUOU 


c:7i n -1- SO 

U ( lU ^ OU 


Sfi -1- 07 

-U.OU ^ U.U/ 


HD 


102870 


72 1-^ 


1400-^ 


OUiZ ^ u^ 


19-1-0 07 
U. -LZ ^ U.U 1 


HD 


121370 


oy .y 


7^C\^ 
1 uu 


cn9S -1- 47"^ 
ouzo ^ ^ 1 


9/1-1- 07 
U.Z^ ^ U.U i 


HD 


139211 


OO.U 


loUU 


«900 -1- SO 
DZUU IC oU 


04 -1- 07 
-U.U^ III U.Ul 


HD 


160691 


yu.u 


zuuu 


UUUU ^ ou 


'?9 -U 07 
U.OZ ^ U.U f 


HD 


165341 


iUi • / 


Af^OOJ 

•yiOuu 


r:QOO -1- SO 


19-1-0 07 
U.IZ III U.U f 


HD 


170987 


t;r: c _|_ n om 
Oo.O ^ u.o 


1 000"^ 

lUUU 


fk^Af) -\- 100"^ 
UU^U ^ lUU 


-u.io ^ u.uu 


HD 


181420 


/ o.u 


lOUU ^ ouu 


UUoU ^ lUU 


00 -1- Ofi° 

u.uu ^ u.uu 


HD 


181906 




1 Q1 9 -L A7P 
±y ±Z ^ '4: / 


UOUU ^ -LUU 


oil -1- 1 4P 

- U . -L 1 ^ U . i 


HD 


186408 


103.1^ 


2150^ 


uozu ^ uu 


n 01 -1- 0*^9 

U.U-L ^ u.uo 


HD 


186427 




9c; crnq 


1 OU iC OU^ 


Of^ -1- 09'3' 

u.uo Hi u.uz^ 


HD 


185395 




9000^^ 
ZUUU 


U 1 ^0 ^ lOU 


Ozl''^ 


HD 


203608 


120 4^ 


9«no^' 
zuuu 


oyyu ^ ou 


74 -1- 07 
-U. i ^ ^ U.U / 


HD 


210302 


RQ '^'^ 

oy .0 




C9Qt; -L QO 
OZOO ^ OU 


01 -1- 07 

U.U-L ^ U.U ( 


KIC 


3632418 


fin ^i'^ -1- n 'iV 
Du.oo m u.o 1 


1110-1- 90 

i±iu in zu 


CI _i_ 70 
DlOU iC (U 


n 1 Q -1- 07 
-u.iy III U.Ul 


KIC 


3656476 


QQ 70 -1- n 99 
yo. / u ^ u.zz 


1 QAf) -1- 9^^ 

ly^iU ^ zo 


r:700 -1- 70 
/ UU ^ / U 


Q9 -1- 07 

U.OZ ^ U.U f 


KIC 


4914923 


QQ fil -1-0 '^9 


1 QQr: _|_ 
±000 ^ UU 


c;Q/in -L 70 
OO'lU ^ / u 


14-1-0 07 

U. -L^ ^ U.U f 


KIC 


5184732 


or: CO _|_ n 


9070 -1- 90 
zu / u ^ zu 


r:Q9c; -L 70 

oozo ^ / u 


'?Q -1- 07 

u.oy ^ U.U ( 


KIC 


5512589 


fiQ CIO -L n 

uo.oz ^ u.oo 


1 9A0 -1- 9^^ 
±z^u zo 


C71 n -L 70 

1 -LU ^ / U 


04 -1- 07 
U.U^ ^ U.U i 


KIC 


6106415 


1 f)'A R9 -1- 9Q 
iuo.Oj^: ^ u.zy 


99Qr: _|_ 90 
zzotj ^ zu 


cQCfi -U 70 
oyou ^ / u 


n 1 1 -L n 07 

-U. 1 -L ^ U.Ul 


KIC 


6116048 


1 no ^A -\- n 99 

-LUU. l-'-t U.ZZ 


91 90 -1- 90 
zizu ^ zu 


CQQc; _|_ 70 
ooyo ^ / u 


n 9f; -L n 07 

-U.ZU ^ U.Ul 


KIC 


6603624 


110 9H -1-0 9'^ 


Z-^rUO Hz OU 


K^iOO -1- 70 

oouu in / u 


9R -1- 07 
U.ZD III U.U 1 


KIC 


6933899 


79 1 -1- n 9^^ 

/ Z. J.tJ ^ U.ZO 


1 Q70 -1- "^0 
±0 1 V ^ OU 


r:Qon 1 70 
OOOU ^ ( u 


n 01 -1- 07 

-U.U-L ^ U.Ul 


KIC 


7680114 


85.13 ± 0.14 


1660 ± 25 


5815 it 70 


0.10 ± 0.07 


KIC 


7976303 


50.95 ± 0.37 


910 ± 25 


6050 it 70 


-0.52 ± 0.07 


KIC 


8006161 


148.21 ± 0.19 


3545 ± 140 


5340 it 70 


0.38 ± 0.07 


KIC 


8228742 


63.15 ± 0.32 


1160 ± 40 


6000 it 70 


-0.15 ± 0.07 


KIC 


8379927 


120.86 ± 0.43 


2880 ± 65 


5960 it 125 


-0.30 ± 0.07 


KIC 


8760414 


116.24 ± 0.56 


2510 ± 95 


5765 it 70 


-1.19 ± 0.07 


KIC 


10018963 


55.99 ± 0.35 


985 ± 10 


6300 it 65 


-0.47 ± 0.50 


KIC 


10516096 


84.15 ± 0.36 


1710 ± 15 


5900 it 70 


-0.10 ± 0.07 


KIC 


10963065 


103.61 ± 0.41 


2160 ± 35 


6015 it 70 


-0.21 ± 0.07 


KIC 


11244118 


71.68 ± 0.16 


1405 ± 20 


5705 it 70 


0.34 ± 0.07 


KIC 


11713510 


69.22 ± 0.20 


1235 ± 15 


5930 it 52 




KIC 


12009504 


88.10 ± 0.42 


1825 ± 20 


6060 it 70 


-0.09 ± 0.07 


KIC 


12258514 


74.75 ± 0.23 


1475 ± 30 


5950 it 70 


0.02 ± 0.07 



Refe rence s: All HD star measurement s without referenc e s are ta ken from | Br untt ct al.l |j2010lV "iTeixeira et aL |j2009h. 'IVauclair et all 
1 I2OO8I). 'iBouchv &: Carried ||2003|). '^Theyenin et al.l 1I2OO5I). teruntt eta l. (2010) -ICarrier et al.l ||2005|). Worth et al.l ||2009^ . 
^Carrier et al.l ll2005h. iBouchv et al.l |j2005|). JjCarrier fc Eggenbergeil 1120061). *lMosser et al.l ||2008|'). iBallot et al.l 1I2OIII ') 



"iMathur et al. 



iBarban et al.l | |2009|) . feruntti l l2009t) . q Garcia et al. I ll2009h . ^Metcalfe et al.l j2012bl '). ICuzik et al.l 1 I2OI j) . 

When (Au) is not explicitly given it is cal culated from the hig hest amplitude / = modes. Values for the KIC 
stars are taken from lMathur et al. 1 1I2OI2I) . 



^Deheuvels et al 



0.025 with a systematic offset of -0.003 ± 0.005, indicating 
no systematic trends. 



5 PRECISION IN logs FOR A LARGE SAMPLE 
OF Kepler STARS OF CLASSES IV/V 

Our primary objective was to test the accuracy of a seismic 
logp by using bright nearby targets that have independent 
mass and radius measurements. We showed in Sect. 13.41 (see 
Fig. [2| that our accuracy should be on a level of 0.01 dex 



with a precision of ~0.02 dex using the set {(Ai/),i/max,7off} 
for the small sample of stars covering the range 3.9 < log g < 
4.6. In this section we investigate the precision in log g for a 
sample of 403 V/IV stars (loggr > 3.75) observed with the 
Kepler spacecraft by employing the same analysis methods. 
In particular we pay attention to systematic errors by 1) 
using different sets of observational constraints, 2) compar- 
ing results using two different methods which incorporate 
different stellar evolution codes and physics, and 3) we also 
show the distribution of errors as a function of magnitude 



10 0. L. Creevey et al. 



Table 5. Derived seismic logg for an extended list of 
Sun-like oscillators, logg and loggno[pe/H] refer to using 
{{A!^),i/niax,T'(,fi,[Fc/H]} and {(A^),!-'max,T'efi} as the constraints 
in the analysis. For a homogenous analysis we adopted 1.1 /zHz 
and 5% as the errors on (A;^) and i/max. We list 2(t uncertainties 
for logg. 



Star Name 



logg 



log3no[Fc/Hl 









(dcx) 






(dcx) 




HD 


10700 


4 


55 


± 





02 


4 


57 


± 





02 


HD 


17051 


4 


40 


± 





02 


4 


39 


± 





03 


HD 


23249 


3 


81 


± 





02 


3 


78 


± 





04 


HD 


49385 


3 


98 


± 





02 


3 


97 


± 





04 


HD 


52265 


4 


28 


± 





02 


4 


24 


± 





02 


HD 


61421 


3 


98 


± 





02 


3 


97 


± 





03 


HD 


63077 


4 


22 


± 





02 


4 


25 


± 





03 


HD 


102870 


4 


11 


± 





02 


4 


10 


± 





04 


HD 


121370 


3 


82 


± 





03 


3 


82 


± 





03 


HD 


139211 


4 


20 


± 





02 


4 


21 


± 





02 


HD 


160691 


4 


22 


± 





02 


4 


21 


± 





02 


HD 


165341 


4 


54 


± 





02 


4 


54 


± 





02 


HD 


170987 


3 


98 


± 





02 


3 


98 


± 





03 


HD 


181420 


4 


15 


± 





02 


4 


15 


± 





03 


HD 


181906 


4 


22 


± 





02 


4 


23 


± 





02 


HD 


186408 


4 


28 


± 





02 


4 


28 


± 





03 


HD 


186427 


4 


36 


± 





02 


4 


35 


± 





03 


HD 


185395 


4 


22 


± 





02 


4 


23 


± 





02 


HD 


203608 


4 


35 


± 





02 


4 


37 


± 





04 


HD 


210302 


4 


23 


± 





02 


4 


24 


± 





03 


KIC 


3632418 


4 


00 


± 





03 


4 


01 


± 





04 


KIC 


3656476 


4 


23 







02 


4 


23 







03 


KIC 


4914923 


4 


21 


± 





02 


4 


21 


± 





03 


KIC 


5184732 


4 


26 


± 





02 


4 


25 


± 





02 


KIC 


5512589 


4 


05 


± 





02 


4 


04 


± 





04 


KIC 


6106415 


4 


29 


± 





02 


4 


30 


± 





03 


KIC 


6116048 


4 


25 


± 





02 


4 


27 


± 





03 


KIC 


6603624 


4 


32 


± 





02 


4 


32 


± 





02 


KIC 


6933899 


4 


09 


± 





02 


4 


09 


± 





04 


KIC 


7680114 


4 


18 


± 





02 


4 


17 


± 





04 


KIC 


7976303 


3 


89 


± 





02 


3 


92 


± 





04 


KIC 


8006161 


4 


48 


± 





02 


4 


48 


± 





02 


KIC 


8228742 


4 


02 


± 





03 


4 


02 


± 





04 


KIC 


8379927 


4 


36 


± 





02 


4 


40 


± 





02 


KIC 


8760414 


4 


31 


± 





02 


4 


35 


± 





03 


KIC 


10018963 


3 


96 


± 





03 


3 


96 


± 





03 


KIC 


10516096 


4 


17 


± 





02 


4 


18 


± 





03 


KIC 


10963065 


4 


28 


± 





02 


4 


29 


± 





04 


KIC 


11244118 


4 


09 


± 





02 


4 


08 


± 





03 


KIC 


11713510 


4 


05 


± 





04 


4 


05 


± 





03 


KIC 


12009504 


4 


21 


± 





02 


4 


21 


± 





03 


KIC 


12258514 


4 


12 


± 





02 


4 


12 


± 





03 



and logp and summarize the uncertainties and systematics 
as a function of magnitude. 



5.1 Observations 

During the first 9 months of the Kepler mission targets 
to be monitored with a 1 minute cadence during 1 month 
each were selected by the KASClfJ. These stars were chosen 



^ Kepler Asteroseismic Science Consortium, see 
|http : //astro . phys ■ an . dk/KASC/ | 
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Figure 4. Comparison of log g derived by Rad ExlO with (top) 
detailed seismic analyses jMathur et al.l l2012l l for 22 Sun-like 
stars observed by Kepler , and (bottom) six common stars from 
ISilva Aguirre et al.l l|2012l 'l. 



based on information availa ble in the Kepler Input Cata- 
log, KIC, (|Brown et al.ll201ll ) and were expected to exhibit 
solar-like oscillations. A total of 588 stars w ith values of 
log q between 3.0 and 4.5 dex were analysed (|Garcfa et al.l 
l201l|) and had their global seismic quantities d etermined 
(|Huber et alll201ll : IVCTner et al.|[201ll : IChaplin et al...201ll ). 
In this paper we concentrate on a subset of 403 less-evolved 
stars with logg values between 3.75 and 4.50 dex derived 
from RadExlO, the range for which we have validated our 
method. 

The global seismic quantities have been determ i ned us - 
ing the SYD pipeline as described by [Huberelai] l|2009t ). 
which uses the reference values of {Av}q = 135.1 /iHz and 
I'maxo = 3,090 /iHz. To avoid systematic errors (Chaplin, 
et al. in prep.), we adopted these same values in our grid. 
The uncertainties on the seismic quantities include a contri- 
bution from the sca tter between different analysis pipelines 
jVerner et al.|[201ll . Chaplin et al. in prep.). Figure [5] shows 
the cumulative distribution for the errors in the seismic ob- 
servations, (Ai/) (top) and Vmax (bottom) in units of /xHz 
and %, respectively, for the 403 stars. Our choice for abso- 
lute and relative errors is for consistency with units used in 
the recent literature. We show the cumulative distribution 
of the errors in order to see the typical errors for 50% and 
80% of the sample, which justifies the errors that we used 
in Sect. |3^ These are less than 1.1/2.0 /xHz for (Av), and 
less than 5%/8% for i^max, respectively. 

The Tcff derived by the KIC have been shown not to 
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Figure 5. Cumulative distributions of the errors in (Av) (top) 
and fmax (bottom) for 403 Kepler s tars with derived log g between 
3.75 and 4.5 dex. Data taken from lHuber et ah I ll201lh . 



be ac curate on a star-to-star basis (|Molenda-Zakowicz et al] 
l201ir). and so the g round-based Kepler s upport photometry 
Brown et al.ll201ll ) was re-analysed by iPinsonneault et al.l 
20121 ) o determine more accurate TcS for the ensemble of 



Kepler stars. These are the temperatures that we adopt for 
our first analysis and we refer to them as Tcffpin. In their 
work they c onsider a mean [Fe/H] = - 0.20 ± 0.30 for all 
of the stars. ISilva Aguirre et all (|2012|) adapt an infra-red 
flux method presented by Casagrande et alj (|2010l ) for de- 
termining stellar temperatures. In their work they apply this 
method to the large ensemble of Kepler stars to provide an 
alternative determination of Tcff. We also adopt their Tcff 
determinations in order to study the effect of biases in tem- 
perature estimates on the derived value of log p. We refer 
to these temperature estimates as Tofiirfm. Support spectro- 
scopic data have also been obtain ed for 93 of the star s and 
the metallicities are presented in iBruntt et al.l (|2012t ). We 
further include these data to study the effects of possible bi- 
ases arising from lack/inclusion of metallicity information. 



5.2 Seismic log 3 from different sets of observations 



Using our validated method described in Sect. 13.21 we cal- 
culated values of logg and their uncertainties for the sam- 
ple of 403 stars using the set of observations comprising 
{T'effpin,(^J^)ii'max} whilc adopting a mean [Fe/H] = - 0.20 
±0.30 dex. We refer to this set as the reference set. The dis- 
tribution of the uncertainties as a function of log g is shown 
in Figure [6] Here it can be seen that typical uncertainties in 
log p for this set of 403 Kepler stars is below 0.02 dex (there 
is one star with an error of 0.05 dex), with a mean value of 
0.015 dex. 

In Figure [7] we show the difference in the fitted log g 
while considering different input observational sets com- 
pared to the reference set 'log grof'. The subsets are; Set 
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Figure 6. Distribution of uncertainties in log g o-(log g) as a func- 
tion of log g for a sample of 403 Kepler stars, using RadExlO with 
observational constraints comprising {T(,ffpin,(Ai/),;/max}. 



1 which considers {/S.v) and Tcftpin only and Set 2 which 
considers (Az/), I'max, and Tefiirfm. We note that for all of 
the analyses [Fe/H] was constrained to -0.20 ± 0.30 dex, 
which corresponds to >90% of the models of the grid. 

Inspecting the top panel of Fig. [7] (Set 1) one can see 
that by omitting i/max as an observed quantity can result 
in differences of over 0.05 dex for a very small percent- 
age of the stars, but the absolute difference between the 
full set of results is -1-0.005 dex with an rms of 0.01 dex. 
Here, we note that several authors have investigated the re- 
lation between (At^) and i/ max and find tigh t correlations 
(|Bedding fc KieldsenI [2OO3I : IStello et al.]|2009l ). The uncer- 
tainties arising from a set of data with less constraints usu- 
ally increases and we indeed find an increase in the uncer- 
tainties (T of ~ 0.01 dex. The extra scatter of 0.01 dex is 
taken care of in the larger a. 

Inspecting the lower panel of Fig. [7] (Set 2) one can see 
that the Tcff derived by using different photometric scales 
results in a mean difference in logg of -0.002 dex (i.e. no 
significant overall effect) with an rms scatter of 0.007 dex. 
This latter fact implies that we can expect to add just under 
0.01 dex to the error budget in logp by considering Tcff 
derived from different methods. We found a similar result in 
Sect. |33] for (3 Hydri. 

The mean value of the derived uncertainties (ct) in log g 
for Set 1 and 2 are 0.023 and 0.015 dex, respectively, while 
those for the reference set are 0.015 dex. The accuracy of 
these log g (if we consider the reference set to be correct) is 
within a precision of Icr. 

Figure |8] compares the derived log g using the refer- 
ence set of observations, to tho se with measured Toff and 
[Fe/H] from lBruntt et all (|2012l ). The absolute mean resid- 
ual is 0.002 dex and is highlighted by the dotted grey line. 
We find that logg can differ by up to 0.02 dex by includ- 
ing [Fe/H]. This 0.02 dex is also consistent with what we 
found in Sect. 13.51 for /3 Hydri when we considered different 
metallicity constraints. 



5.3 Comparison of results using different codes 

To investigate the possible source of systematics arising 
from using a different evolution code and input physics, 
we determined \ogg using a second pipeline code, Yale- 
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Figure 9. The difference between tiie vaiues of logg obtained 
by RadExlO and Yale-Biiam as a function of RadexlO log g and 
normalised by the Yale-Bham errors. See Sect. [531 for details. 



Figure 7. Differences in log g using subsets of observational con- 
straints. The reference set comprise ((Ai/), i/max, T^sviri)- Set 1 
and 2 comprise ((A^), Tcffpin) and ((Av), fmax, TeffiR), respec- 
tively. For all sets [Fe/H] = -0.20 dex. The itlcr error bars are 
average error bars measured over bins of 0.1 dex corresponding 
to Set X. See Sect.[52]for details. 
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Figure 8. Differences in the derived seismic logg from RadExlO 
between adopting [Fe/H] = -0.20 ± 0.30 for all of th e stars (log g) 
and a dopting [Fe/H] from spectroscopic analyses bv lBruntt et al.l 
||2012|) . logg[po/H]- See Sect. [5l2]for details. 



Bham (jGai et al.ll2"oTl] : ICreevev et al.ll2012al l. Details of the 
code can be found in the cited papers. Here, it suffices to 
know that the method is very similar to t hat of RadExlO, 
but t he evolution code is based on YREC (|Demaraue et al] 
I2OO8I ) in its non-rotating confi guration, with the foUow - 



ing specifications: OPAL EOS feogers fc Navfonovl I2OO2I') 
and OPAL high-temperature opacities (jlglesias fc Rogers! 
1996t) supplemented w ith low-temperature opacitites from 



Ferguson et al. (|2005|). a nd the NACRE nuclear reaction 
rates ( Angulo et al.lll999l V Diffusion of helium and heavy- 
elements were included, unlike RadExlO. 

Figure [9] shows the difference in the derived 
value of logg between RadExlO and Yale-Bham using 
{(Ai^),i^niax,T'offPin,[Fe/H]=-0.2}. Wc show the normalised 
difference, i.e. divided by the Yale-Bham errors, which are 
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Figure 10. Distribution of derived uncertainties in the seismic 
logg from RadExlO for 403 Kepler V/IV stars as a function of 
SDDS r magnitude. 



very similar to the RadExlO errors. As can be seen from the 
figure the agreement in loggr between the different methods 
is within Icr or < 0.01 dex for 98% of the stars and 2a for 
all stars except one. A small absolute difference of 0.005 is 
found with an rms of 0.005 (units of a). This is most likely 
due to the different physics adopted in the codes. A very 
slight systematic trend is present with a slope of -0.005 ± 
0.005. The agreement between the results using the different 
codes and methods is very encouraging. We can consider an 
upper limit of 0.01 dex as a typical systematic error arising 
from employing different codes. 



5.4 Uncertainties and systematic errors in the 
seismic \ogg as a function of magnitude 

In Figure [TU] we show the uncertainties in logp as a function 
of SDDS r magnitude (obtained from the KIC) while using 
the reference set of observations. As the figure shows, preci- 
sion improves with apparent brightness, where more reliable 
measurements have smaller error bars. We can expect typ- 
ical uncertainties in loggr of less than 0.03 dex for the full 
sample with a mean value of < 0.02 dex. 

In Table |6] we summarize for different apparent mag- 
nitude bins 1) the mean uncertainties (a) and 2) the mean 
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Table 6. Summary of uncertainties a and systematic errors s in 
the seismic log g from the Kepler sample of stars as a function of 
magnitude. We give the mean values () of each over the sample 
that falls into the corresponding magnitude bin. The systematic 
errors s have subscripts which refer to possible sources of errors 
in the input T^jf, [Fe/H], and the code that is used to generate 
the stellar models (see Sects. [5T2l and [5731 for details). In the final 
row we summarize the mean uncertainties and adopted systematic 
errors for the full sample. 



r 






(«[Fc/H]) 


(•^codc) 


(mag) 


(dex) 


(dex) 


(dex) 


(dex) 


6 < r <9 


0.012 


0.004 


0.014 


0.005 


9 < r < 10 


0.014 


0.004 


0.009 


0.004 


10 < r < 11 


0.016 


0.009 




0.006 


11 < r < 12 


0.017 


0.009 




0.006 


6 < r < 12 


0.015 


0.01 


0.02 


0.01 



systematic errors (s) arising from adopting different tem- 
perature scales (sTcff), metallicities (s[po/H]) and diflerent 
pipeline codes (scodo). Note that the sample size for the bin 
with r < 9 is much smaller than the others and this is the 
reason we find a larger S[p(,/Hi for r < 9. In the last row 
we give the mean uncertainty and the adopted systematic 
errors over the full sample. 

5.5 Summary of results for the sample of Kepler 
stars 

From the sample of Kepler data with 3.75 < logg < 4.50 
derived from seismic data, we obtain a mean precision o of 
between 0.01 and 0.02 dex (max 0.03 dex for 99% of stars) 
when (Ai/), i^max, and Tcff are used as the input constraints 
for the pipeline analysis. A typical systematic error a of 0.01 
dex can be added to account for a possible systematic error 
in TeS, such as that obtained by applying different calibra- 
tion methods to photometry. Similarly, we can add a sys- 
tematic error of 0.02 dex due to an incorrect or lack of a 
[Fe/H] measurement. We also showed that using different 
models and physics in the pipeline analyses yields results 
in logg consistent to within 0.01 dex i.e. almost no model 
dependence. 

Eliminating one of the seismic indices from the set of 
data yields a zero offset in the results for the full sample of 
stars and a typical scatter of 0.01 dex, although differences 
of over 0.05 dex were found for < 1% of the stars. 



6 DISCUSSION 

The first objective of this study was to investigate if we can 
determine logg reliably using global seismic quantities and 
atmospheric data. We showed that our method is reliable by 
comparing our results with values of logg derived from di- 
rect mass and radius estimates of seven bright nearby stars. 
We then applied our method to a list of 40 Sun-like stars 
and derived logg to wit hin la (~ 0.01 dex ) of th ose from 
iMathur et al] (|2012l ) and lSilva Aguirre et al.l (|2012l ) . We also 
showed that for a sample of 400-f Kepler stars (6 < r mag < 



12) with logg between 3.75 and 4.50 dex typical uncertain- 
ties (J of less than 0.02 dex can be expected and we estimated 
systematic errors s of no more than 0.04 dex arising from 
errors in TcS and [Fe/H] measurements and using different 
codes. 

All of the Kepler stars will be observed by the Gaia 
mission, and for this reason they provide a valid set of cali- 
bration stars, by constraining log g with precisions and accu- 
racies much better than spectroscopic or isochrone methods 
provide for the current list of calibration stars for Gaia. The 
astrophysical parameters to be determined from Gaia data 
using the astrometry, photometry and BP/RP spectropho- 
tometry are Tcff, Aq, [Fe/H], and logg. By ensuring that the 
GSP_Phot methods deliver logg consistent with the seismic 
value will reduce the uncertainties and inaccurcies in the 
other parameters. Moreover, an independent logg provides 
an extra constraint for the determination of the atmospheric 
parameters from GSP_Spec where degeneracies between Toff 
and logg severely inhibit the precision of atmospherically 
extracted parameters and chemical abundances. 

While in this paper we concentrated primarily on us- 
ing Kepler data to determine logg, we note that both the 
CoRoT and Kepler f ields have gr e at pot ential in other as- 
pects. For example, iMiglio et al.l (|2012l ) and iMiglio et al.l 
(|2013l ) determine distances, masses, and ages of red gi- 
ants from the CoRoT fields to constrain galactic evolution 
models, while ISilva Aguirre et al.l (|2012h develop a method 
which couples asteroseismic analyses with the infrared flux 
method to determine stellar parameters including effective 
temperatures and bolometric fluxes (giving angular diame- 
ters), and hence di stanc es. The distances ob t ained by both 
IMiglio et ai] (|2013l ) and ISilva Aguirre et all (|2012l ) can be 
compared directly with a Gaia distance for either calibra- 
tion of Gaia data or stellar models. By combining data from 
CoRoT and Kepler with those from Gaia we can also deter- 
mine extremely precise angular diameters by coupling a seis- 
mically determined radius with the Gaia parallax. Finally, 
these stars will also be excellent calibrators for the FLAME 
workpackage of Gaia, which aims to determine masses and 
ages for one billion stars. 



7 CONCLUSION 

In this work we explored the use of F, G, K IV/V stars ob- 
served in the Kepler field as a possible source of calibration 
stars for fundamental stellar parameters from the Gaia mis- 
sion. Our first objective was to test the reliability /accuracy 
of a seismically determined log g, and using a sample of seven 
bright nearby targets we proved that seismic methods are ac- 
curate (see Fig. [2| by obtaining results to within 0.01 dex 
of the currently accepted logg values. We showed, however, 
that the accuracy of the input atmospheric parameters does 
play a role in the accuracy of the derived parameters. For /3 
Hydri we found that errors in the atmospheric parameters, 
and in particular [Fe/H], can change logg by 0.02 dex. 

We then applied our grid-based method RadExlO to 
an extended sample of stars from the literature. We de- 
rived their seismic logg and these are given in Table [5] 
We showed that a grid-based logg is consistent with the 
values of log q obta ined by detailed seismic analysis from 
iMathur et aiT(|2012h for 22 stars and in excellent agreement 
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with the results from ISilva Aguirre et al.l (|2012l ) for the 6 
common stars. We find typical precisions in \ogg of 0.02 
dex. 

We finally studied the distribution of errors in log g from 
two analysis methods for a sample of 403 Kepler stars with 
3.75 < logg < 4.5, and we obtained a typical uncertainty of 
between 0.01 and 0.02 dex in \ogg for F, G, K IV/V stars 
(see Table [B]). We can add a total of 0.04 dex as a system- 
atic error which arises from the adopted temperature scales 
(0.01 dex) and measured metallicities (0.02 dex) as well as 
the grids of models used (0.01 dex), which differ in evolu- 
tion codes and input physics. The precisions in the data are 
unprecedented, and the systematic errors are much smaller 
than those stemming from any other method, esp ecially for 
7 < V < 12 stars (see e.g. ICreevev et al.|[2012af Fig 1 and 
Table 3 which compare spectroscopically derived logg for 
five Kepler stars with V ~ 11). 

As a final remark, we highlight that while there are 
>500 IV/V stars in the Kepler Held and some in the CoRoT 
fields that exhibit Sun-like oscillations, there are also lOOO's 
of red giants in both CoRoT and Kepler fields with these 
same measured quantities, however, the accuracy of these 
stars' seismic \ogg is yet to be studied. 
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